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@) Arithmetic operation apparatus for elementary function. 



An arithmetic operation apparatus performs an arithmetic operation of a drfferentiable elementary 
function including a sine function sin(x) v a cosine function cos(x), arc tangent function arctan(x), an 
exponential function e* a logarithm function log # (x), an Inverted number 1/x, a square root Vx, and an 
inverted number of a square root 1/Vx. The apparatus comprises a divider receiving an initial value (x) 
for dividing the initial value into a more significant digit portion H and a less significant digit portion L (x 
= H + L), and 2^word memories of "k" banks receiving the more significant digit portion H of the initial 
value (x) and storing a value obtained by multiplying the more significant digit portion H by a previously 
calculated constant b 0 = f(H) or b k = (1/k f) x f« (H) (where k is posffive integer). An arithmetic circuit 
composed of a multiplier and an adder receives the less significant digit portion L of the initial value (x) 
and an output of the memories for executing the following polynomial : 

f(x) = bo + L x { Dl + L x • . . (bfc., + b k x L} • • •) 
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ARITHMETIC OPERATION APPARATUS FOR ELEMENTARY FUNCTION 



Background of the Invention 
Field of the invention 

5 

The present invention relates to an arithmetic operation apparatus for elementary functions in computers.. 
Description of related art 

10 Computers for handling numerical data has to have numerical functions, particularly elementary functions 
including a sine function sin(x), a cosine function cos(x). arc tangent function arctan(x), an exponential function 
e* a logarithm function log,(x), an inverted number 1Ac, and a square root Vx. As a manner for calculating these 
elementary functions, the following algorithms. have been known : >r : 

(1) approximate value method 

15 This method utilizes a polynomial approximation such as phebyshev expansion, Taylor expansion, and 
others, and a rational formula approximation such a rorttirfu^jtection expansion, which have been ordinarily 
used for obtaining an elementary function. However, which of the expansions should be selected dependentty 
upon a function to be obtained. In addition; in the' ca$e of c§ldulation of the elementary function of a floating 
point number in a doubled; precision (precision of 52 digits Wfi bjnajy number), multiplication and addition must 

20 be repeated more than ten times, and therefore, a longlarithmetic operation time ts required. 

(2) CORDIC method - \ j \ ' \ ) . ^ \' \ 

\ This Is used for obtaining a trigonometric function and an inverted trigonometric function. This can unitarily 
calculate a plurality of elementary functions. However, for calculation of the elementary function In a doubled 
precision (precision of 52 digits in a binary number), addition/subtraction and shift must be executed 52 x 3 
25 times, and therefore, a long arithmetic operation time is required; , ' 

(3) STL method , V. > . v ■ _ } 'c 

J ' This Is used foran exponential function and a logarithm function. However, for calculation of the elementary 
function in a doubled precision (precision of 52 digits in a binary number), addition/subtraction and shift must 
"be executed 52 x 2 times; and therefore, a long arithmetic operation time is required. 
so (4) Newton method' * 7 ; ~ jj 

f ; This can be used only for a square root and a cube root, and a calculation equation should be selected 
dependency upon a function to be obtained. An initial value is given in the form of an approximate value, and 
adaltion/subtractiori/rnultiplication/division is repeated until a requited degree of precision is obtained. If th ini- 
tial value has; a relatively high precision, the repetition can be completed one or two times. For calculation of 
35 the elementary function in a doubled precision (precision of 52 digits in a binary number), the method can be 
realized with one time of reading of a table ROM (read only memory) and three to six times of addition/sub- 
traction and multiplication. Therefore, a required arithmetic operation time Is relatively short 

: In any of ithe above mentioned conventional methods, it is rare that the same algorithm can be applied to 
a full range of argument Therefore, it is an, ordinary practice to convert the argument so that calculation is xec- 
40 uted in a limited range. ^ j ; / . - ; - \ , \ 

For example, in the case of calcutation of sin 8, it Is assumed that 9 = p x n/4 + x (where p is integer), and 
6in(x) and cos(x) are calculated in a range of 0 ^ x < */4. A procedure of limiting a range of argument as in this 
method Is called an "argument reduction". 

As seen from the above description, the above mentioned elementary function arithmetic operation 
45 methods are disadvantageous in the following points : 
(1 ) The arithmetic operation time is long. 

Multiplication and addition must be repeated more than tep times in the approximate value method, and 
shift and addition/subtraction must be repeated more than ten times In the CORDIC method and the STL 
method. Therefore, the arithmetic operation time is very long. * 
so (2) A plurality of elementary functions cannot unitarOy be hand(ed. 

An optimum approximate. value Is different if the elementary function Is different, and therefore, the 
approximate valu method cannot unitarily process a plurality of elementary functions. On the other hand, the 
STL method and the Newton method are limited in elementary: functions which can be treated. The CORDIC 
method can unitar&y handle the trigon metric function and the inverted trigonometric function, but cannot treat 
55 the exponential function and th logarithm function. 
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Summary of the Invention 

Accordingly, it is an object of the present invention to provide an elementary function arithmetic operation 
which has overcome the above mentioned defect of the conventional one. 

Another object of the present invention is to provide an elementary function arithmetic operation having a 
table for an elementary function f(x) and previously calculated numbers fW (x)/k I which are obtained by multi- 
plying a k-order differentiation of f(H) with a predetermined constant, the elementary function arithmetic oper- 
ation capable of quickly and precisely obtaining the elementary function arithmetic operation f(x) by a one-tim 
searching of the table, "k" times of multiplication and *k" times of addition (where k is an integer in a rang of 
1 to 6 inclusive) 

The above and other objects of the present invention are achieved in accordance with the present invention 
by an arithmetic operation apparatus for arithmetically operating a differentiate elementary function including 
at least one of a sine function sin(x), a cosine function cos(x), arc tangent function arctan(x), an exponential 
function e* a logarithm function 3og e (?:), an inverted number a square root ^c, and an inverted number of 
a square root 1A^x, comprising dividing means receiving an initial value (x) for dividing the initial value into a 
more significant digit portion H and a less significant digit portion L (x «= H + L), 2«-word memories of "k* banks 
storing previously calculated constants ^ = f(H) or b* = (1/k !) x fW(H) (where k is positive integer) which are 
obtained by multiplying a k-order differentiation of f(H) by given constants, each of the memories receiving as 
an address the more significant digit portion H cf the initial value (x) for outputting one of stored previously cal- 
culated constants designated by fths address, and at least one arithmetic circuit composed of a multiplier and 
an adder and receiving the less significant digit portion L of the initial value (x) and an output of the memories 
for executing the following polynomial : 

f(x) = bo + L x {bf.+'Lx - - - (bfc.t + bfcx L) * • •) ■, 

In one embodiment, assuring k = 3, th© elementary function f(x) can be obtained with precision of 4n digits 
by three times of multiplication and three times of addition. In another embodiment, assuming k *= 3, the elemen- 
tary function f(x) can be obtained with preccscon of 2n digits by one time of multiplication end one time of addition. 

Now, a principle of the etementary function arithmetic operation in accordance with the present inv ntion 
wfll be described. 

Assume that the elementary function f(x) is infinitely differentiate Li a rang© of a given argument Th 
elementary ft/notion f(x) can he developed into a Taylor series. Assuming x - a, the Taylor development can 
be expressed as follows : 

f(a+5)*f(a) +f (l) (a)x5/l! 

+ f (2) (a)x8/2! 
+ f< 3 >(a) x 8/3 ! 



+ fW( a )x8/k! (i) 

This Taylor development is disadvantageous in that if an absolute value of a is large, an error becomes 
large, and the order number "to? of the approximation equation becomes Barge. On the other hand, if the absolute 
value of 6 is made very small, the practalon of the arithmetic operation can be elevated, and the order number 
"k" of the approximation equation crn be reduced to an range of 1 to 5. 

The elementary function arithmetic operation apparatus in accordance with the present invention is Inten- 
ded to make the range of 5 as narrow as possible and to fall the order number to a possible extent, by pre- 
viously calculating values erf f(a), fW(a), C^(a), ... f*>(a) foras many iniUai numbors "a" as possible In th equation 
d). 

Here, consider th case in which an elementary function i{x) is sought in a binary number of 4n digits, from 
a given number V represented by a binary number of 4n digits (0 s x< 1). 

First, the given number V is divided into a more significant digit portion H and a less significant digit portion 
L by units of "n" digits. H and L correspond to "a" and "5" of the equation (1), respectively. 
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4n 

x = Z {x k x 2' k } « H + L (2) 
where xk = (1, 0} 

0 < X < 4 ^-iA^ ^-y t )^i.r:. (3) 



H = Z i( x k x2-M . ^ . (4) 



L^X [ xkKfc^ . (5)' 



Hera, putting a = H and 6 = L in the equation (1). since L < 2* and Li 4 < 2 4n , the term of L 4 and subsequent 
terms can be neglected or ignored; Therefore, the equation (1) can be obtained with sufficient precision by limi- 
ting as k = 3. Namely, rO.c-'e-.r'j f^-. : -:'.: . 
SO fW^H^KH^L*^^ r ■ ^ 

If this equation (6) is optimized In order to minimize the' number of multiplication operations, the following 
equation (7) can be obtained : 

f<x) = 1(H)*UIF 1 W r^r • 

In this equation (7) T if there is prepared a table ROM which receives H as an input and outputs the following: 



bo = f(H) 

b k = f (k >(H)/k! (8) 



(where k = 1 f 2, 3) 
the foflowing equation (9) can be obtained : 

f(x) = bo+Lx^t+LxCba+Lxbj)} (9) 
Therefore, f(x) can be obtained with three multiplication operations and three addition operations. 
45 Assuming that f(x) has no singular point within a designated range of V, the value of b k will never become 

large. In addition, since L<2«. L*<2-*\ and L*<2-*\ the constant table for the equation (9) is sufficient if it has 
4n digits for bo, 3n digits for b 1( 2n digits for bj, and n digits for b 3 . 

Therefore, in order to obtain a value of the elementary function f(x) with precision of 4n digits, only a table 
ROM of (4n+3n+2n+n) digits x 2" words, multipliers and adders are required. Here, since the precision of mul- 
so tiplication of bs*L can be limited to W digits; soma of the required multipliers and adders can have a reduced 
precision. In other words, all the required multipliers and adders are not required to have the same degre of 
precision of 4n digits, . • 1 u~ ? iwi''" * -v.: 

in addition, some elementary function f(x) has an argument range of 1^x2. However, since an integer 
portion of V Is fixed to "1\ the table ROM Is sufficients it has 2 n words. 
55 The above and other objects, features and advantages of the present invention will be apparent from the 
following description of preferred embodiments of the* invention with ref rence to th accompanying drawings. 

Brief Description of the Drawings 

Figure 1 Is a block diagram of an embodiment of the elementary function arithmetic operation apparatus 
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in accordance with the present invention ; and 

Figure 2 is a bl ck diagram of another embodiment of the elementary function arithmetic operation 
apparatus in accordance with the present invention. 

Description of the Preferred embodiments 

Referring to Figure 1 , there te c ;wn a block diagram of an embodiment of the elementary function arithme- 
tic operation epparatus in accordance with the present invention. The shown apparatus includes a register 1 
for receiving and holding an input variable "x". This register 1 has a capacity of 4n bit or more. More significant 
"n" bits of a decimal fraction portion of the variable V held in the register 1 are read out as a H signal 2, and 
least significant "3n" bfSs of tfie variable V held in the register 1 are road ott as a L signal 3. More significant 
•2n" bits of th© Seast significant "Sn" bits are outputted as a 12 signal 4, and more significant "n" bits of the least 
significant "3n* bits are outputted as a L1 signal 5. 

In addition, tho shown apparatus includes a ROM (read only memory) 6 of 4n bits x 2" words for holding 
values of f(H) t a ROM 7 of 3n bite x 2" words for holding values of f *>(H), © ROM 8 of 2n b!ts x 2° words for 
holding values of ftf>(H) / 2, and a ROM 0 of n bits x 2" words for holding values of *3>(H) / 6. Furthermore, the 
shown apparatus includes a multiplier 10 of "n" bits x V bits receiving the L1 signal 5 and an output of the 
ROM 9 for the purpose of calculating a product of the L1 signal 5 and the output of the ROM 9 t and an adder 

11 of 2n bits for adding an output of the ROM 8 with a data which is obtained by rightward shifting an output 
of the multiplier 10 by "n* bits. An output of the adder 11 and the L2 signal 4 are supplied to another multiplier 

12 of "2n" bits x "2n" bits, which in turn outputs a product of the output of the adder 1 1 and the L2 signal 4. An 
output of the multiplier 12 is rightward shifted by "n" bits, and inputted to one input of another adder 13 of 3n 
bits, which has the other input connected to receive an output of the ROM 7. An outpift of the adder 13 and the 
L signal 3 are supplied to a third multiplier 14 of "3n" bits x "3n" bits, which in turn outputs a product of the 
output of the adder 13 and the L signal 3. An output of the multiplier 14 is rightward shifted by V bits, and 
inputted to one input of a third adder 1 5 of 4n bits, which has the other input connected to receive an output f 
the ROM 6, so that the n-blt rightward-shifted output of the multiplier 14 and the output of the ROM 6 are added. 
An output of the adder 15 is connected to an output register 16 of 4n bits. 

The following table shows various elementary functions f(x) which can calculated by the apparatus shown 
in Figure 1 and which therefore are to b© held in the ROM 6, coefficients to be held in the ROMs 7, 8 and 9, 
and the range of the argument V. 
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Table 1 



J&lwUIvliUUY 

Function 
f(x) 


ISl - KJTOCT 

Differentiation 
CO 


2nd - ureter 
DefferentiaUQ Ql 
f<2)(x)/2, 


3rd - Order 
Differentiation/3 ! 
> r,, f<3) (x)/6 


Range of 
Argument 


sin (x) 


cos (x) 


-sin (xy2vrf:> 


I * :- COS (30/6 ' 

i .-■■■..»!;:" -'; ";;'->j| Mi, 


0<x£Tt/4 


cos(x) 


— sin(x) -.. 


-*cos (x)/2 ' f ^ J 


: ; ■ sin <x)/6 

\ i' t ' i ;■'■».*..■ t .(.' 




arctan (x) 








0£x<1 


; <« 2 + J) 




3 (x 2 - l) 3 , 


c x 


■ y 




cV6 


0£x£Log e 2 


Iog„ (x) 


l/x 


1/2* 






l/x 


-l/x* 


l/x 3 


-l/x 4 


l<x<2 


fx" 


1 

2/F 


-1 

8xVx" 


1 

16x 2 /x" 


l<x<4 


l//x~ 


-1 
2x VT 


3 

8x 2 /3T 


-5 
16x 3 /x* 


l£x£4 



For example, if the elementary function f(x) to be obtained is sin(x), it could be understood from the table 
1 that sin(H) is stored in the ROM 6 f cos(H) is stored In the ROM 7 and — ein{Hy2 and — cos(H)/6 are stored 
in the ROMs 8 and 9, respectively. 

In the apparatus shown in Figure 1 t if the argument V is inputted and registered in the input register 1, 
the H signal 2. which is composed of the more significant V bits of the decimal fraction portion of the argument 
V held in the register 1, is supplied to the ROMs 6 to 9 as an address. On the other hand, the multiplier 10 
executes multiplication of the L1 signal 5 and a valu of ff»>(H)/6 read out of the ROM 9. A value of LxW^/e 
outputted from th multiplier 10 is rightward shifted by V bits for decimal point matching, and is added with a 
value of f&(H)72 road from th ROM 8 f by means of the adder 11. The multiplier 12 multiplies an output value 
of the adder 1 1 by the 12 signal 4. An output of the multiplier 1 2 is rightward shifted by p rt" bits for decimal point 
matching, and is added with a valu of f«>(H) read from the ROM 7, by means of the adder 13. The multiplier 
14 multiplies an output value of the adder 13 by th L signal 3. An output f the multiplier 14 is rightward shifted 
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by "n" bits for decimal point matching, and is add d with a value of f(H) read from the ROM 6, by means of the 
adder 15. 

5 As seen from the above, the arithmetic operation of the elementary function f(x) requires only a total times 

of an access time for the table ROMs and a time required for the three multiplication operations and the thre 
addition operations. 

!n the case cf obtaining a value of a given elementary function f(x) with precision of 52 bits (n=1 3), 

ROM 6 has a memory capacity of 52 bits x 8192 words; * 
10 ROM 7 has a memory capacity of 39 bits x 8192 words ; 

ROM 8 has a memory capacity of 26 bits x 8192 words ; 

ROM 9 hes © memory capacity of 13 bits x 8192 words ; 

Multiplier 10 is a multiplier, cf 13 bits x 13 bits ; 

Multiple 12 is a multiplied 26 bits x 26 bits ; , 
is Multiplier 14 3s a multiplier of 39 bits x 39 bits; . 

Adder 11 is an adder of 26 bits ; * . v - 

Adder 13 is an adder of 39 bits ; and 

Adder 15 is am adder of 52 biSs 0 

In the above case, a total memory capacity of the ROMs is 1,064,960 bits, which can be realized in the 
20 form of an LSI according to a recent semiconductor technique. 

Furthermore, Sf each of the multipliers 10, 12 and 14 is constituted of a carry save adder and a carry propa- 
gate adder, the amount of hardware and the time of arithmetic operation can be effectively reduced by con- 
stituting sach of the multipliers 10, 12 and 14 and the adders 11, 13 and 15 by the cany save adder and by 
putting cany propagate adder after th© adder 15^ % >.-•; 
25 Turning to Rguro 2, there is shown a block diagram of another embodiment of the elementary function 

arithmetic operation apparatus In accordance with the present invention. The first embodiment has a high deg- 
ree of precision for the purpos® of a doubled precision floating point data (such as 52 digits in a binary notation). 
However the second embodiment is for a single precision floating point data (such as 24 digits in a binary nota- 
tion). Specifically, the first embodiment has the precision of 4n digits under k=3 in the equation (3), and on the 
30 other hand, ftho second embcdimant has the precision of 2n digits under k=1. 

In the second embodiment therefore, the equation (9) can be expressed as follows : ,: 

f(x) = bb+Lxb 1 (10) , 
Therefore, f(x) can be obtained with one multiplication operation and one addition operation. In addition, 
sine© H<1 and L<2+ 9 the constant table for the equation (10) is sufficient if it has 2n digits for and "n" digits 
35 fbrb,. • ■ . - • , .. . .. ....... ,, 

The second embodiment includes a register 1a of 2n bits for receiving and holding an input variable V. 
More significant "n" bits of a decimal fraction portion of the variable V held in the register 1a are read out as 
a H signal. 2, and Coast srgnfficant V bits of the variable V held in the register 1 are read out as a L signal 3a. 
In addition, the shown apparatus includes a ROM 6a of 2n bits x 2" words for holding values of f(H), and a ROM 

40 7a of n bits x 2" words for holding values of fM(H), which is a first-order differentiation of f(H). Furthermore, the 
shown apparatus includes a multiplier 10 of "n* bits x "n" bits receiving the L signal 3a and an output of the 
ROM 7a for the purpose of calculating a product of the L signal 3a and the output of the ROM 7a. and an adder 
11 of 2n bits for adding an output cf the ROM 6a with a data which is obtained by rightwarti shifting an output 
of the multiplier 10 by "n* bits. output of the adder Sa cs connected to an output register 16a of 2n bits. 

45 If the argument V is inputted and registered in the input register 1a, the H signal 2. whSch is composed of 

the more significant *n" bits of the decimal fraction portion of the argument V held in the register 1a. is supplied 
to the ROMs 6a and 7a as an address. On the other hand, the multiplier 10 executes multiplication of the less 
significant W bits of the argument V (the L signal 3a) and a value of fW(H) read out of the ROM 7a. A value 
of Lxfl D(H) outputted from the multiplier 10 is rightward shifted by "n" bits for decimal point matching, and is 

so added with a value of f(H) read from the ROM 6a, by means of the adder 1 1 . As the result, the adder 1 1 outputs 
f(x)*f(H) + Lxf(i)(H). 

As seen from the abov , the arithmetic operation of the elementary function f(x) requires nly a total times 
of an access tim for the table ROMs and a time required for the n multiplication operation and the one addi- 
tion operation. 

65 In the case of btainingavalue fa given elementary function f(x) with precision cf 24 bits (n=12), the ROM 

6ahasamemorycapacityof24bitsx4096words ;th ROM 7a has a memory capacity of 12 bits x 4096 words; 
the multiplier 10 is a multiplier of 12 bits x 12 b*s ; and the adder 11 is an adder of 24 bits. 

Therefore, a total memory capacity of the ROMs Is 147,456 bits. Therefore, In order to realize the fght 
elementary functions shown in the table 1, a memory capacity cf 1,179.648 bits is required, which can be 
realized in the form of an LSI. 
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As seen from the above, the elementary function arithmetic operation apparatus in accordance with the 
present invention is advantageous in the following two points : 
5 (1 ) The arithmetic operation time Es short 

In the first embodiment of the elementary function arithmetic operation apparatus, the required processing 
time is a total time of the reading time of the table ROM, a triple of the multiplication time and a triple f the 
addition time. In the second embodiment, the required processing time is a total time of the wading time of the 
table ROM, the multiplication time and the addition time. 
10 Now, assuming that the reading time of the table ROM is 0.20 us, the multiplication time is 0.20 fxs, and 
the addition time is 0.05 |is, the required processing time is 0.9S us in the first embodiment and 0.45 [is in the 
second embodiment :.jyp>:j:].i,:\f-. i....-- 

(2) A plurality of elementary functions f(x) can be unitarily handled: 

The elementary function arithmetic operation apparaius in accordance with the present invention can 
is change over the elementary function f(x) by changing data stored in the table ROM. Therefore, if table ROMs 
corresponding to the eight kinds of elementary functions shown iri the Table 1, one of the eight kinds of elem n- 
tary functions, sin(x) f cos(x), erctan(x), e", log»(x), 1/x. Vx. and i/sfcc ban be selectively calculated. 

The invention has thus been shown and described with reference to the specific embodiments. However, 
t should be noted that the present invention is in no way limited tc the details of the illustrated structures but 
20 changes and modffications may be made within the scope of the appended claims. 

■ Claims - :] [ * - f - " 

se ■ ?1, An arithmetic operation apparatus fo^ 

ing dividing means receiving an initial value (x) for dividing said initial value into a more significant digit - i 
pcrtimi H and a iess significant digit portion L (x « H + L) f 2"-wQrd memories of "V banks storing previously f 
caicidated constants bo s f(H) er « <i,ic !) x fW<H> (where k is positive Integer) which are obtained by ^ 
muitipiyirig a k-order differentiation of f(H) by given constants, each of said memories receiving as an # 

so address said mora significant digit portion H of said Initial value (x) for outputting one of stored previously 

calculated constants designated by said address, and an arithmetic circuit composed of at least one mul- 
tiplier and at least one adder and receiving said less significant digit portion L of said initial value (x) and • ft 
an output of said memories for executing the following polynomial : - J " * ■ n 

f(x) = ho + L x {bi ♦ L x . - • (b^i + b|c x L) • •) .< 

36 ■ -'■ ■ '" • ■ - ' v : * " : i)s ' : ' [/] \ : * : ' 

Z, An aiithrrtetic operation apparatus claimed in Claim 1 wherein said k is 3 and saW arithmetic circuit is com- ^ 5 
posed of three multipliers and three adders so that a value of the elementary function can be obtained with 

precision of4n digits. < - ,!; - -' ■ :;r: '' r f\ ^•■ V! - ■■' 

40 3c An arithmetic operation apparatus claimed in Claim 1 wherein said k is 3 end said arithmetic circuit is com- 
posed of one multiplier and one adder 30 that a value of the elementary function can be obtained with pre- 
■ ciston of 2n digits. -■■ : : - - : -- r -'' 

4; An arithmetic operation apparatus claimed in Claim 1 wharein said elementary function includes at least 
45 one of sine function sin(x) f a cosine funcSon cos(x), arc tangent function arctan(x), an exponential function 

e*, a logarithm function log c (x), an inverted number 1/x; a square root Vx. and an inverted number of a 
square root 1/>/x. < >* ... ■ « % : . U*.a ; . *. ./ \ 
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FIGURE 2 
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